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Abstract: The impact of the system design on the control of coupling 
between planar waveguide modes and surface plasmon polaritons (SPP) is 
analyzed. We examine how the efficiency of the coupling can be enhanced 
by an appropriate dimensioning of a multi-layer device structure without 
using additional gratings. We demonstrate that by proper design the length 
of the device can be dramatically reduced through fabrication a surface 
plasmon resonance sensor based on the SPP-photon transformation rather 
then on SPP dissipation. 
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1. Introduction 

Recent advances in the application of surface plasmon polaritons in the fields of chemical and 
bio-sensors and nano-science has stimulated booming research and development. Surface plas- 
mon resonance (SPR) based optical bio-sensors exploit a special type of the electromagnetic 
surface waves coupled to the collective electron excitation — surface plasmon polaritons (see 
e.g. [1-6]) to probe interactions between an analyte under study and a bio-molecular recogni- 
tion element immobilized on the SPR sensor surface. One of the most important characteristics 
of the SPPs is their ability to detect small scale effects and field interactions at the interfaces 
between a metal and a non-metal. The surface plasmon polaritons attract a great deal of atten- 
tion as a quite unique possibility of the electromagnetic field localization and corresponding 
substantial enhancement of the electric field near the surface. Such a high level of a field lo- 



calization results into the enhanced sensitivity of the SPR sensors. Specifically, SPPs in sand- 
wiched structures have attracted much interest as a promising way to implement surface plas- 
mon resonance optical sensors [7-11]; waveguide-coupled SPR sensor was first demonstrated 
by Harrie Kreuwel in [12]. Many of the key properties of such sensors, using multi-layer sand- 
wiched structures, are determined by the interaction between a waveguide and SPP modes. The 
progress in the area is very fast and some products are already commercially available. There 
is broad range of possible schemes and configurations, and despite successful first demonstra- 
tions of such devices, many design problems are still open and in this paper we examine some 
of these issues. 

We would like to point out that SPP-based technology has a very broad spectrum of potential 
applications. A feasibility of exploiting the high spectral dispersion achievable in plasmonic 
devices in optical communication has been recently discussed in [13]. Surface waves play an 
important role in many fundamental resonant phenomena, such as e.g. the "Wood anomalies" 
in the reflectivity and transmissivity [1, 14, 15] of periodically corrugated metal samples, and 
are already practically exploited in a wide range of practical devices. A very important area 
of the SPP applications is in miniaturization and integration of optical circuits. A possibility 
of high localization of electromagnetic fields using SPP and subwavelength metal structures 
makes it possible to overcome the diffraction limit of light waves which is a critical issue in 
the development of photonic circuits. Considerable efforts have been recently placed into the 
development of plasmon subwavelength guiding and coupling devices [16-19]. 

A promising direction in the development of SPP-based devices and techniques is to com- 
bine surface plasmon approaches with well developed dielectric waveguide technologies, in 
particular with fiber-optics [20-23]. The conventional modern SPR fiber sensors utilize Bragg 
gratings or long period gratings to couple light to the SPPs. Recently the coupling efficiency 
between the core mode of the waveguide and SPP modes of the metal film in planar and cylin- 
drical geometry for sensing purposes have been examined in [24-27]. A variety of possible 
configurations of multi-layered structures offers a range of potential device designs and this 
field is still far from been fully explored. Efficiency of the coupling between SPP modes and the 
waveguide is the key to unlocking the full potential of such sandwiched multi-layer structured 
devices. In this work we analyze the basic properties of the SPP and waveguide mode coupling 
and consider a possibility to optimize the SPR-sensors by system design through tuning device 
parameters rather than through additional gratings written in the waveguide. Understanding of 
the details of the interaction of the SPP and waveguide modes is a critical point in the de- 
velopment of the high sensitivity real-time SPR- sensor and our goal here is to analyze such 
interaction for rather general configuration that can be relevant to many practical situations. 

2. Fields representation and boundary conditions 

Consider eigenmodes of the multilayer structure depicted in Fig. [T] that is typical for the SPP- 
based devices. The indices c, b, w, — and + in what follows relate to the corresponding layers 
as shown in Fig. [T] Taking into account that SPP is a TM-polarized wave and, therefore, in the 
geometry under consideration it can interact only with TM-polarized waves, we consider TM- 
polarized eigenmodes (E = {E x ,0,E z }, H = {0,if,0}) of the structure. The time dependence 
is moved out of the equations by the corresponding Fourier transform and analysis of the plane 
waves having the form exp(— icot). We are interested here in solutions that present outgoing 
or decaying waves both in the upper half-plane z<0(t=— ), and in the region z > L + i + d 
(t = +) for any given frequency CO and tangential component of the wavevector k x , 

H T (x,z)=H T exp{ik x x + iTk T z [z-S T ^(d + £ + L)}}, (1) 



where k — co/c is the vacuum wavenumber, k\ = \/k 2 £ T — k%, Im(fcJ),Re(fcJ) > 0. Similarly, 




Fig. 1. Geometry of the problem and a typical dielectric permittivity distribution. The struc- 
ture of the SPP modes for both thick and thin metal films are schematically shown along 
with a waveguide mode. 



H M (x,z) = £ H M \° cxp{ik x x + iok¥[z-z M ]}, 

G=± 



(2) 



in the metal (M — c, z G [0, J]), the buffer (M = b, z G + and the waveguide (M = w, 
z G [J + d + £ + L]) layer. Here z c = 0, = d, z w — d + £ and the branch of the square roots, 
is chosen so that Re(^) > 0, lm(k^) > and thus the sign a = +(— ) corresponds to the 
waves propagating or decaying in the positive (negative) direction of the z-axis, respectively. 
The x-dependence in all the layers is the same due to conservation of the tangential momentum 
(or, equivalently, as imposed by the boundary conditions). The electric field components can 
be derived from the Maxwell equations. Assuming all the media shown in Fig. [T] to be non- 
magnetic we have 

E?(x,z) = -^-^-H m (x,z), (3) 



ke m dz 



(4) 



Jc£jfi dx 

where m = T,c,b,w, and £ m denotes dielectric permittivity in the m-th layer (medium), and 

El(x,z) =E*exp{ik x x + iTkl[z-8 Ti +(d + £ + L)}}, (5) 

E?(*,z)= L E^ G ^{ik x x + iGk^[z-z M ]}. 



(6) 



G=± 



Here u = x, y stays for the corresponding components of the electric field. The amplitudes of 
the electric field, , can be expressed by means of Eqs. ([3]), ^ and ([T]), ^ as 
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Introducing dimensionless variables, 



a = k x /k, P" 



z 

k£ m 



(9) 



Eqs. ([7]), ^ can be presented in the form 



a 



(10) 



E M\a = Gj3 M H M\a E M\o^_^ R M\a^ 

Continuity of the tangential components of the field crossing the interfaces results in the sys- 
tem of eight linear homogeneous algebraic equations for the eight amplitudes of the magnetic 
field, H M \ G , H T . The nontrivial solutions of these system (eigenmodes of the layered struc- 
ture) exist when the determinant is equal to zero. This leads after straightforward algebra to the 
dispersion relation: 



where 



J8 w [ia w tan - b w ] + j8 + [a w - ib w tan 0] = 0, 



ob 

a w = a b + b b tanhF, b w = £- \a b tanhF + b b ] , 



(12) 



a b = l + ^-tanh<D, b b = ^ 
% = jS e is the surface impedance of the metal film, and 



<& = —ik c d, <j> 



KM 



tanh<£- 



(13) 



(14) 



Resolving Eq. fl2| ) yields the dispersion curves, CO = co(k x ) [or k x = K(co)] corresponding to 
different modes of the system. For the sake of clarity, but without loss of generality, in what 
follows we neglect frequency dependence of all dielectric permittivities except that one of the 
metal layer, £ c = £ c (co). For the latter we will use the Drude model, £ c (co) = £oo — (6^/co((0 + 
iy), where is high-frequency value of the permittivity £ c (g>), co p is the plasma frequency, 
and 7 is the damping rate. In Table [T] the typical parameters for the commonly used materials 
(for design of the sensors) are shown, here n m = y/£^, m = —,b,w,+, denote refractive indices 
of the corresponding dielectric layers. 



Table 1 . Table of common material parameters for the SPR sensing devices* 
metal silver [7] gold [7, 24, 25] aluminium [28j 
dielectric core /i = 1.585 [29] n= 1.476 [28] n= 1.47 [24] 
dielectric cladding n = 1 .439 [29] n = 1 .473 [28] n = 1 .45 [24] 
* For optical constants of silver, gold and aluminium see [30-32]. 



The field amplitudes can be expressed in terms of the upper medium magnetic field ampli- 
tude, H~\ 

\ P J z (15) 

H w\o = ^_ | fl w _ ab y^ H + = H - [ fl w cos _ ib « sin ft . 

The dispersion relation includes a number of limiting cases corresponding to different com- 
binations of thick or thin limits of the building layers. For instance, in the limit of a large 



thickness of the buffer layer, i — > oo (F ^> 1), the dispersion relation splits into two dispersion 
branches, 

P b + P+ ~ i (V + ^) tan* = 0, (16) 

d* + &* = 0, (17) 

which correspond to the waveguide TM modes and SPP modes of the metal film, respectively. 
In turn, the derived SPP dispersion relation, Eq. ( [17] ), for a rather thick metal film (<£ ^> 1) limit 
leads to standard SPP modes, a b + b b = => (P b + £ ) (j3 " + £ ) = 0. In the case of a finite value 
of <I> the dispersion relation describes the hybridized double-interface-localized SPP modes, 

(j3^tanh<I> + ^)(j3-tanh<I> + ^)-^ 2 cosh- 2 <I> = 0. (18) 

More specifically, when the dielectric permittivities of the superstate and the buffer coincide, 
£b = £_, then Eq. ( [T8] ) splits into two branches and yields the dispersion relations for the short- 
range (symmetric) SPP and long-range (antisymmetric) SPP modes [1], see Fig.[T] 

pb = pb^ _£ tanh (<D/2), j5 b = l5 b = -§coth(*/2), (19) 

where subscripts "+" and "— " correspond to the long-range (LR), and short-range (SR) mode, 
respectively. 

Similar to that, the solitary waveguide modes (£ — > oo) depend on the waveguide layer thick- 
ness and permittivity. Namely, in this limiting case the solution of Eq. ( fT6| ) reads 

+ p> n = 0,±l,... (20) 

At the crossing points of the dispersion curves the phase matching conditions are satisfied 
and an interaction between waves of different branches is most efficient. However, in order to 
stress the difference between the situations when the phase matching conditions are provided 
by an additional grating and the problem considered here we will use in what follows term 
hybridization of modes. In the problem considered here the finite thickness of the buffer layer 
results in hybridization (combination) of the waveguide and SPP modes. The resulting hybrid 
eigenmodes of the system might have spatial distributions with energy both in the dielectric 
waveguide and in the SPP modes. We would like to emphasize that the most interesting situa- 
tions for the design of the sensors takes place when the level of the mode hybridization (strong 
coupling of the waveguide and SPP modes) is rather high. This occurs when the buffer layer is 
relatively thin and in the vicinity of the points in the co — k x plane where the dispersion curves 
of the free SPP and waveguide modes intersect. In the vicinity of such points the hybridization 
is high and control of operation near this point can be provided by the adjustment of the buffer 
layer. This consideration outlines the design rules and explains the choice of the parameters in 
the examples presented below. 

The solutions of the dispersion relation of the multi-layer system, Eq. ([12]), are plotted in 
Fig. [2] where we use the normalized frequency parameter kL = Lco/c, versus the squared nor- 
malized tangential component of the wavenumber (i.e., the squared effective refraction index, 
n 2 ff = a 2 = c 2 k 2 /(Q 2 ). The most interesting area for the parameters lies is the vicinity of the 
intersection between the silver- superstate SPP branch and that of the lowest-frequency (funda- 
mental) waveguide mode, see inset in Fig. [2] 

The coupling of these modes due to the finite-thickness of the buffer layer results in a repul- 
sion of the dispersion branches. However, in the vicinity of the inoculating curves intersection 



£ w kL = —— arctan 




Fig. 2. Dispersion relations for the multi-layer system of Fig.[T]Thick (red) curves represent 
the solutions of the complete dispersion relation given in Eq7p2]>. Thin (gray) WG-named 
curves correspond to the dispersion relation of the solitary dielectric waveguide, Eq. {20|, 
and a thick silver metal film, d^°°. The dashed (blue) curves represent the high-frequency 
and low-frequency solutions of the dispersion relation for a solitary thin silver metal film. 
The parameters are: n- = 1.37, n^—n^ — 1.439, n w = 1.585, the thicknesses of the layers 
are d = 58 nm, £ = 400 nm, L = 250 nm. 



the distance between the modified dispersion branches is relatively small. The intersection de- 
fines the "resonance frequency", k res L = 3.37 (fijfes/(2rc) = 6.44 • 10 14 Hz) for L = 250 nm, see 
Fig. [2] that optimizes the considered structure properties for sensor applications. Note that our 
approach can be also used in the reverse problem - when the design of the multi-layer structure 
is adjusted to some required operational frequency. 

The magnetic field evolution along the x-axis for the solitary metal film SPP modes, i.e., 
when the buffer and the waveguide thicknesses are zero- valued, i,L — > 0, could be derived 
from Eqs. ([I]), ^ and ( fl"5] ), where k x and CO obey Eq. fl8] ). This distribution in plane is 
shown in Fig. [3 for the above "resonance frequency", COr QS . 

Field evolution for the eigenmodes of the solitary dielectric waveguide (we consider here the 
limit or vanishing buffer and the metal layer thicknesses: i — > 0) could be derived from the 
same equations but for k x and CO obeying Eq. p0| ). The corresponding magnetic field distribu- 
tion is shown in Fig.|4| for CO = (Ores. 

The magnetic field distributions in the plane (z,x) for three eigenmodes of the multilayer 
system at the "resonance frequency" are shown in Figs. [5] [6] and [7] 

It is seen that by an appropriate fitting of the parameters of a multi-layer structure it is pos- 
sible to switch between different regimes and to control efficiency of coupling between the 
waveguide mode and SPP modes. 

3. Transmission properties of the multi-layer structure 

Consider now transmittance of the described above multi-layer system with the length of the 
metal layer that differs from the dielectric waveguide length. We calculate here the transmit- 
tance as a ratio between the total output energy flux at x = 100 /im, and the input energy flux 
at x = 0. The input energy flux at the starting point x — corresponds to the solitary waveguide 
fundamental mode, and the output energy flux is obtained by integrating the x-component of 
the energy flux density across the section x = 100 jim from z = — 3 fim to z = 6 jim. The 
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Fig. 3. The magnetic field distribution for the two SPP eigenmodes of the solitary metal film 
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as in Fig. [2] The a (b) figure is the field distribution for the high-frequency (low-frequency) 
SPP mode of the metal film corresponding to the upper (lower) SPP dispersion curve in 
Fig.0 
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Fig. 4. The magnetic field evolution for the solitary waveguide mode, a 
parameters are the same as in Fig. [2] 
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electromagnetic field within the structure was obtained by direct numerical simulations of the 
problem by means of the finite element method (FEM). The calculation length of 100 jim 
was chosen because it was found to be enough to observe the effect. The computation of the 
longer structures requires significantly larger computer resources. However, it does not add sig- 
nificantly new information — use of a longer computational domain leads to a more precise 
characterization of a fine structure of the region near the minimum at the resonance frequency 
(it becomes slightly wider and deeper), but the main features of the frequency dependence of 
the transmittance are practically not changed. 
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Fig. 5. The magnetic field evolution for the eigenmode of the system corresponding to 
the point 1 in Fig. [2] Compare the field distribution with that of the right in Fig. [5] The 
parameters are the same as in Fig. [2] 
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Fig. 6. The magnetic field evolution for the eigenmodes of the system corresponding to the 
points 2 (a) and 3 (b) in Fig. [2] 



The results of the numerical modeling of the transmission coefficient for the long metal film 
length > 100 fim (formally, — > oo) are presented in Fig. [8] The frequency range in the 
calculations corresponds to the stripe area in the Fig. [2] 

One can see rather narrow dip in the transmission coefficient near the defined above "res- 
onance frequency", shown by the vertical line. Evidently, the transmission falls due to exci- 
tation of the SPP modes accompanied by strong enhancement of the dissipative losses. Other 
local minima of the transmission coefficient are caused by appearance of the radiation losses, 
cf. Fig. [2] We would like to emphasize that in the considered structure there exists 3 eigenmodes 
of the system simultaneously (for some given input frequency), see Fig. [2] This fact leads to the 
field interference pattern easily seen in Fig. [9] For the "resonance frequency" the largest period 




Fig. 7. Squared magnetic field evolution for three eigenmodes of the system for the points 
1, 2 and 3 of Fig. [2] and the dielectric permittivity distribution. The parameters are the same 
as in Fig. [2] 



of the pattern is given by the following simple expression: 

2n 

^beat 7 j — , (21) 

where k x 2 and k x 3 are the real parts of the propagation constants of the long-period modes 2 and 
3 of the complete system, and could be found from the dispersion relation. In this particular 
instance Xb ea t ~ 69.5 jim. The small features in the pattern are due to interference of the short- 
period mode 1 with the long-period modes 2 and 3. That is, corresponding characteristic length 
of the pattern is X^ eat ~ 2n/(k x \ —k X 2) — 3.9 /im, cf. Fig. [9] These features are in excellent 
agreement with the results of the recent experimental work [27], where authors consider similar 
geometry and discuss the coupling length and energy exchange between waveguide and SPP 
modes. 

At the input of the system main part of the energy was in the dielectric waveguide and then 
is gradually transferred to the metal film. It is seen that the beat half-period Xb ea t/2 gives us a 
characteristic length for energy exchange between the core- and the high-frequency metal film- 
localized mode (coupling length) and at the distance x = X^ at /2 the maximal transfer of the 
energy into the metal-film-localized mode occurs, see Fig. [9] Consequently, low transmission 
coefficients (corresponding to almost total energy transfer from waveguide to the SPP) can 
be achieved by choosing the metal film length of the order J2f opt = Xb ea t/2, see Fig. [8] This 
distance is defined as a minimal length of the metal film that can provide for the total energy 
transfer from the waveguide core to the metal. Since the process of the energy transfer from 
the core to the SPP and in the backward direction in the absence of losses is cyclic, there 
exist a number of metal film lengths for which the energy total vanishing in the core appears: 
j£f m = Xbeat/2 + fft^beat; m = 0,l,2,.... However, the minimal metal film length has advantage 
of a large contrast between resonant and non-resonant transmission (transmission at resonant 
and non-resonant frequencies) due to the small dissipation losses in the shorter metal film. 
Additionally, the SPP propagation length in our geometry at the resonance frequency is about 
100 |Um, that is longer than the minimal length of the metal film. In other words we propose 
here to make a SPR sensing device, that is shorter (comparable) than SPP propagation length, 
in contrast to the commonly studied devices with metal film that is much longer than SPP 
propagation length [10, 1 1 ,36] . We anticipate that devices with such a minimal length might also 
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Fig. 8. Frequency dependence of the transmission coefficient for the metal film length > 
100 jUm (solid line), and for j£f = 40 jUm ~ X^ QSLt /2 (dashed line) obtained by direct numeric 
simulations. The ripples are caused by the finite film length. The arrows a and b indicate the 
frequency values below which the radiative losses arise for the modes 3 and 2, respectively. 
Dotted lines present results obtained by the mode expansion procedure. The thickness of 
the buffer layer is £ = 400 nm. 




Fig. 9. The squared magnetic field distribution at the resonance frequency, direct numeric 
simulation. 



have enhanced sensitivity. In the recent work [37] it was observed that with decrease in sensing 
region length, the sensitivity, signal-to-noise ratio are both increasing and the resolution also 
improves. However, analysis of the sensitivity and signal-to-noise ratio issues for the devices 
with minimal sensing region length proposed here are beyond the scope of this paper and will 
be discussed elsewhere. 

We stress that it is clearly seen in Fig.[8]that the effect is even more pronounced as compared 
with the long-film system. Indeed, the transmittance is about zero at the resonance frequency 
for the dashed curve corresponding to the metal film length of j£f « This means, in 

particular, that there is no need to make a millimeter-long [10, 11, 36] film in order to observe 



strong transmittance suppression due to the SPR resonance. This feature might be of a particular 
interest for device miniaturization and integration. The ripples in the frequency dependence 
of the transmittance coefficient are due to the finite length of the metal film. In Fig. [10| the 
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Fig. 10. The squared magnetic field evolution at the resonance frequency for the structure 
with the metal film, which is cut at some special length of 40 /im. Direct numeric simula- 
tion. 



distribution of the squared magnetic field is depicted. It is simulated at the resonance frequency 
for the system with the metal film of the length j£f = 40 /im. It is seen that, indeed, near the 
metal film end the energy flux is concentrated close to the film, demonstrating approximately 
a total energy transfer from the core into the SPP mode and consecutive scattering out of the 
system at the metal film edge. As a result, the energy flux in the cross-section x = 100 /im is 
practically disintegrated. 

Finally, we would like to discuss briefly the possibility to apply approximate approach that 
can be used when a direct numerical modeling becomes time consuming. The fields evolu- 
tion and transmittance through the structure can be approximately found using the well-known 
mode expansion method [35]. We would like to point out that there exist rather advanced the- 
oretical methods for analysis of the excitation of surface plasma waves in integrated-optical 
waveguide structures [36], however, in our geometry it appears to be sufficient to apply a more 
simple and straightforward mode expansion approach in the following realization. Note that 
in the considered configurations the homogeneity of the multi-layer structure along the x-axis 
is preserved (e.g. no gratings are used), there will be no energy transfer (energy exchange) 
between the eigenmodes of the complete structure during the propagation. Therefore, it is pos- 
sible to express any arbitrary incoming field distribution — the fundamental waveguide mode 
of the dielectric waveguide H w {x,z, co) in our example — at the point x — in the basis of the 
eigenmodes of the complete system, Hi(x,z, (o), see Fig. [7] In other words, we calculate what 
energy each mode would carry through the structure. Applying this expansion we obtain the 
field evolution in an arbitrary point x along the structure: 

H(x,z,(D) = ^Ci((D)-Hi{x,Z,(0), 

r k-(co) (22) 

Q(«) = /tf w (0^ 

where Hi are the eigenmodes of the complete system, H w is the input fundamental waveguide 
mode of the dielectric waveguide, k x i are the wavevectors of the eigenmodes of the complete 
system. Using Eq. ( [22] ) it is easy to calculate the transmittance through the structure, Fig. [8] 



The comparison of the direct numeric simulations and the modes expansion method shows an 
excellent agreement as it can be seen in Fig. [8] 

4. Conclusion 

We have examined the properties of a multi-layer hybrid device combining the waveguide and 
surface plasmon polariton modes in a planar geometry. Our main focus here was on the analysis 
of the feasibility to control mode coupling efficiency through the system design rather than via 
using additionally written gratings. We have demonstrated that by the proper design of the 
multi-layer structures it is feasible to achieve substantial reduction in the length of the device. 

We have shown that in the considered scheme (without using additional gratings) of the 
SPR-based sensors there exists an optimal length of the metal film. Such an optimal length is 
defined by the characteristic (coupling) distance of the energy exchange between the waveguide 
core and SPP mode. Moreover, there is a set of preferable metal film lengths (sensing region 
lengths), which are multiples of the coupling distance. The minimal sensing region length offers 
advantage of higher sensitivity and better signal-to-noise ratio and an improved resolution as it 
was observed in the recent work [37]. Minimization of the sensing region length also provides 
a better contrast between resonant and non-resonant transmission (transmission at resonant and 
non-resonant frequency), due to the small dissipation losses in the shorter metal film. In general, 
SPR sensor with the minimal metal film length utilizes the effect of SPP-photon transformation 
rather then SPP dissipation in the metal plate. In other terms our design analysis provide a 
guidance where to cut a metal film to make a short sensing device. The obtained results might 
be important for miniaturization and integration of a range of nano-optic systems such as e.g. 
bio- and chemical sensors, frequency- selective polarizers, and other devices. 

We have shown that an appropriate choice of the operational (resonance) frequency could 
be derived from a rather simple examination of the intersections of the waveguide and SPP 
dispersion curves resulting in their hybridization. This opens a way of an efficient excitation of 
the surface plasmon resonance by the waveguide mode and vice versa. In the situations when 
the required operational frequency (wavelength) is defined by the external conditions our re- 
sults can be used as the design rules for a specific multi-layer device to operate at the desirable 
frequency. Evidently, hybridization — coupling of the waveguide and SPP modes occurs not 
only at a single frequency point, but rather in some spectral interval. Shift of an operational fre- 
quency from an optimal value leads to gradual reduction of the efficiency. Our results indicate 
that a characteristic length of the devices satisfying the resonance conditions (accompanied by 
a high level of hybridization) required for an efficient energy transfer between waveguide and 
SPP modes could be as short as tens of micrometers. This is quite an important departure from 
the characteristic millimeter-scale [10, 11,36] lengths of conventional devices. Our scheme can 
be potentially generalized to other geometries, for instance, it can be applied to the cylindrical 
(fiber) geometry. By an appropriate selection of the parameters the operational (resonance) fre- 
quency could be shifted into the near infrared region including telecom windows near 1500 nm. 
This can be achieved, for instance, by increasing the waveguide layer thickness up to 350 nm 
and simultaneous decrease of the metal layer thickness down to 30 nm. The minimum transmis- 
sion of the multi-layer structure considered here depends essentially on the dielectric properties 
of the media adjacent to the metal film (the superstate). Therefore, small variations in the re- 
fractive index of this medium would result in a strong shift of the transmission dip offering new 
opportunities and perspectives for the sensor applications. 
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